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. We discuss here a self-consistent method to calculate the properties of the cold asymmetric nuclear 

^Nj ' matter. In this model, the nuclear matter is dressed with s-wave pion pairs and the nucleon- 

nucleon (N-N) interaction is mediated by these pion pairs, u> and p mesons. The parameters of 
^ ', these interactions are calculated self-consistently to obtain the saturation properties like equilibrium 

binding energy, pressure, compressibility and symmetry energy. The computed equation of state is 
then used in the Tolman- Oppenheimer-Volkoff (TOV) equation to study the mass and radius of a 
neutron star in the pure neutron matter limit. 
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PACS numbers: 21.65.+f,21.30.Fe, 24.10.Cn,26.60.+c 



J;^ ! I. INTRODUCTION 



The search for an appropriate nuclear equation of state has been an area of considerable research interest because 
of its wide and far reaching relevance in heavy ion collision experiments and nuclear astrophysics. In particular, the 
studies in two obvious limits, namely, the symmetric nuclear matter (SNM) and the pure neutron matter (PNM) have 
'■ helped constrain several properties of nuclear matter such as binding energyper nucleon, compressibility modulus, 
T^j- [ symmetry energy and its density dependence at nuclear saturation density po to varying degrees of success. Of 

late, the availability of flow data from heavy ion collision experiments and phenomenological data from observation of 
^ compact stars have renewed the efforts to further constrain these properties and to explore their density and isospin 
. content (asymmetry) variation behaviours [1, @, 0| • 

■ One of the fundamental concerns in the construction of nuclear equation of state is the parametrization of the 
OO , nucleon-nucleon (N-N) interaction. Different approaches have been developed to address this problem. These methods 
• can be broadly classified into three general types [8|, namely, the ab initio methods, the effective field theory approaches 
^ . and calculations based on phenomenological density functionals. The ab initio methods include the Brueckner- 
• '-j ■ Hartree-Fock (BHF) d, M, El approach, the (relativistic)Dirac-Brueckner-Hartree-Fock (DBHF) 0, HI, 0, [11, ^ 
' calculations, the Green Function Monte-Carlo (GFMC) pvl ITsl method using the basic N-N interactions given 
^ by boson exchange potentials. The other approach of this type, also known as the variational approach, is pioneered 



by the Argonne Group [20|, [2lj . This method is also based on basic two-body (N-N) interactions in a non-relativistic 
formalism with relativistic effects introduced successively at later stages. The effective field theory (EFT) approaches 
are based on density functional theories [22, 23] like chiral perturbation theory [24, 25]. These calculations involve a 
few density dependent model parameters evaluated iteratively. The third type of approach, namely, the calculations 
based on phenomenological density functionals include models with effective density dependent interactions such as 
Gogny or Skyrme forces [1^ and the relativistic mean field (RMF) models [13, [H, 0, [s^] • The parameters of these 
models are evaluated by carefully fitting the bulk properties of nuclear matter and properties of closed shell nuclei to 
experimental values. Our work presented here belongs to this class of approaches in the non-relativistic approximation. 

The RMF models represent the N-N interactions through the coupling of nucleons with isoscalar scalar a mesons, 
isoscalar vector lo mesons, isovector vector p mesons and the photon quanta besides the self- and cross-interactions 
among these mesons [1^ . Nuclear equations of state have also been constructed using the quark meson coupling model 
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(QMC) [3l| where baryons are described as systems of non-overlapping MIT bags which interact through effective 
scalar and vector mean fields, very much in the same way as in the RMF model. The QMC model has also been 
applied to study the asymmetric nuclear matter at finite temperature [s^l ■ 

It has been shown earlier [sl, [s^l , that the medium and long range attraction effect simultcd by the a mesons in RMF 
theory can also be produced by the s-wave pion pairs. This "dressing" of nucleons by pion pairs has also been applied 
to study the properties of deuteron[35l| and ^He [sB]- On this basis, we start with a nonrelativistic Hamiltonian density 
with TriV interaction. The repulsion and the isospin asymmetry part of the NN interaction are parametrized by 
two additional terms representing the coupling of nucleons with the lo and the p mesons respectively. The parameters 
of these interactions are then evaluated self-consistently by using the saturation properties like binding energy per 
nucleon, pressure, compressibility and the symmetry energy. The equation of state (EOS) of asymmetric nuclear 
matter is subsequently calculated and compared with the results of other independent approaches available in current 
literature. The EOS of pure neutron matter is also used to calculate the mass and radius of a neutron star. We 
organize the paper as follows: In Section II, we present the theoretical formalism of the asymmetric nuclear matter 
as outlined above. The results are presented and discussed in Section III. Finally, in the last section the concluding 
remarks are drawn indicating the future outlook of the model. 



II. FORMALISM 

We start with the effective pion nucleon Hamiltonian 

7^(x) = 7^^(x) + H„^(x) + Hm(x), (1) 
where the free nucleon part T-In{x) is given by 

HAr(x) = V'^(x) V(x), (2) 

the free meson part 'Hm{'x.) is defined as 

Hm(x) = l[^l + (V^O • (V^O + mV-] , (3) 
and the ttN interaction ^33| is provided by 

iG 



V^(x). (4) 



In equations ([2]) and ([4]), represents the non-relativistic two component spin- isospin quartet nucleon field. The 
single particle nucleon energy operator is given by = {A'P - V^)^/2 with nucleo n mass M and the pion-nucleon 
coupling constant G. The isospin triplet pion fields of mass m are represented by if. 

We expand the pion field operator ipi (x) in terms of the creation and annihilation operators of off-mass shell pions 
satisfying equal time algebra as 

'/'^(x) = ^^(aj(x)''' +a,(x)), .^,(x) = iW^^(a,(x)t - aj(x)), (5) 

with energy uj^ = {m? — V^)^/'^ in the perturbative basis. We continue to use the perturbative basis, but note that 
since we take an arbitrary number of pions in the unitary transformation U in equation ([7]) as given later, the results 
would be nonperturbative. The expectation value of the first term of TLmti'^) ^l- O vanishes and the pion pair 
of the second term provides the isoscalar scalar interaction of nucleons thereby simulating the effects of a-mesons. A 
pion-pair creation operator given as 

= ly /(k) a,(k)t a,(-k)tdk, (6) 

is then constructed in momentum space with the ansatz function /(k) to be determined later. 
We then define the unitary transformation U as 

L/^ef^'-^), (7) 
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and note that [/, operating on vacuum, creates an arbitrarily large number of scalar isospin singlet pairs of pions. 
The "pion dressing" of nuclear matter is then introduced through the state 

I/) = [/|mc) =e(^'-^)|mc), (8) 

where U constitutes a Bogoliubov transformtion given by 

C/t a,{k)U = (cosh /(k)) a,(k) + (sinh /(k)) a,(-k)t, (9) 

We then proceed to calculate the energy expectation values. We consider N nucleons occupying a spherical volume 
of radius R such that the density p — N/{^ttR^) remains constant as {N, R) oo and we ignore the surface effects. 
We describe the system with a density operator /5jv such that its matrix elements are given by (ssj 

p„/3(x, y) = Tr[pN ^/3(y) Va(x)], (10) 

and 

Tr[pNN] = J p,„(x, x)dx ^N^pV. (11) 
We obtain the free nucleon energy density 

(/|^r[p^.H^.(x)]|/> = V 2^ I A/+ ---^ 1 . (12) 




In the above equation, the spin degeneracy factor j — 2, the index r runs over the isospin degrees of freedom n and 
p and kj represents the Fermi momenta of the nucleons. For asymmetric nuclear matter, we define the neutron and 
proton densities p„ and pp respectively over the same spherical volume such that the nucleon density p = pn + Pp- 
The Fermi momenta are related to neutron and proton densities by the relation /cj — (Gn'^pr/j)^ ■ We also define 
the asymmetry parameter y = (pn — Pp)/p- It can be easily seen that pr = |(1 ± y) for t = n,p respectively. 

Using the operator expansion of equation ([5]) , the free pion part of the Hamiltonian as given in equation ^ can be 
written as 

7Ym(x) = aj(x)^ Lo^ ai(x). (13) 
The free pion kinetic energy density is given by 

hk = (/|7^m(x)|/> - ^ j dk ^(k) sinh^ /(k), (14) 



where ^(k) — Vk^ + rn^ ■ Using ex — M in the nonrelativistic limit, the interaction energy density hint can be written 
from equation (g]) as 

h,nt - {f\TT[pM nU^)]\f) ^ ^(/l : ¥'.(x)^.(x) : |/). (15) 
Using the equations ([7]), ([8]) and ([9]), we have from equation (fT5|) 



C^p f 2, f dk /sinh2/(k) . ...mi 
h^nt = ^ 7^ / -TTT +smh2 /(k) . (16) 



2M V(27r)3 J w(k) V 

The pion field dependent energy density terms add up to give hm{— hk + hint) which is to be optimized with respect 
to the ansatz function /(k) for its evaluation. However, this ansatz function yields a divergent value for h^. This 
happens because we have taken the pions to be point like and have assumed that they can approach as near each other 
as they like, which is physically inaccurate. Therefore, we introduce a phenomenological repulsion energy between 
the pions of a pair given by 



3a 
(2^ 



(sinh-' /(k)) e^'-''^ dk, (17) 
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where the two parameters a and Rt^ correspond to the strength and length scale, repectively, of the repulsion and 
are to be determined self-consistently later. Thus the pion field dependent term of the total energy density becomes 
hm = hk + hint + ^m- Then the optimization of with respect to /(k) yields 

r^n 1 

tanh2/(k) = ^ . (18) 

The expectation value of the pion field dependent parts of the total Hamiltonian density of eqn. ([1]) alongwith the 
modification introduced by the phenomenological term becomes 



. with the integrals (t = n,p) given by 



Ir = 



PrJn + Ppip (19) 



(20) 



and w — ti-'(k). 

We now introduce the energy of w repulsion by the simple form 

K = A^p2^ (21) 

where the parameter \^ corresponds to the strength of the interaction at constant density and is to be evaluated 
later. We note that equation (|2T|) can arise from a Hamiltonian density given in terms of a local potential vr{'k) as 

Hfl(x) ^(x)tV'(x) J vr{^ - y)^(y) V(y)rfy, (22) 
where, when density is constant, we in fact have 



The isospin dependent interaction is mediated by the isovector vector p mesons. We represent the contribution due 
to this interaction, in a manner similar to the w-meson energy, by the term 

hp - Xppl (23) 

where p^ = [pn — pp) and the strength parameter Apis to be determined as described later. Thus wc finally write 
down the binding energy per nucleon Eb of the cold asymmetric nuclear matter: 

Eb = --M, (24) 
P 

where e = (^/ + hm + h^ + hp) is the energy density. The expression for e contains the four model parameters a, 
i?^, and Ap as introduced above. These parameters are then determined self-consistently through the saturation 
properties of nuclear matter. The pressure P, compressibility modulus K and the symmetry energy Esym are given 
by the standard relations: 

P - (25) 
dp 

K ^ (26) 



2 9y2 



y=0 



The effective mass M* is given by M* = Af + Vs with Vs = [h^ + h^)/p. 
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III. RESULTS AND DISCUSSION 



We now discuss the results obtained in our calculations and compare with those available in literature. The four 
parameters of the model are fixed by self-consistently solving eqs. ([M)) through ^7} for the respective properties of 
nuclear matter at saturation density po = 0.15 fm~^. While pressure P vanishes at saturation density for symmetric 
nuclear matter (SNM), the values of binding energy per nucleon and symmetry energy are chosen to be —16 MeV 
and 31 MeV respectively. In the numerical calculations, we have used the nucleon mass M = 940 MeV, the meson 
masses m = 140 MeV, = 783 MeV and rup = 770 MeV and the tt - TV couphng constant C^/Att = 14.6. In 
order to ascertain the dependence of compressibility modulus on the parameter values, we vary the K value over a 
range 210 MeV to 280 MeV for the symmetric nuclear matter {y = 0) and evaluate the parameters. It may be noted 
that this is the range of the compressibility value which is under discussion in the current literature. For K values in 
the range 210 MeV to 250 MeV, the program does not converge. The solutions begin to converge for compressibility 
modulus K around 258 MeV. We choose the value 260 MeV for our calculations. In Table [J we present the four 
free parameters of the model for ready reference. 



TABLE I: Parameters of the model obtained by solving the equations (|24p - (|27p self consistently at saturation density. 



a i?^ Ap 
(MeV) (fm) (fm^) (fm^) 
16.98 1.42 3.10 0.65 




FIG. 1: The binding energy per nucleon Eb as a function of relative nucleon density p/ po calculated for different values of the 
asymmetry parameter y. The values y — 0.0 and 1.0 correspond to symmetric nuclear matter (SNM) and pure neutron matter 
(PNM) respectively. 

For this set of parameter values the effective mass of nucleons at saturation density is found to be M* /M = 0.81. In 
the Fig. [TJ we present the binding energy per nucleon Eb calculated for different values of the asymmetry parameter 
y as a function of the relative nuclear density p/po- The values y = 0.0 and 1.0 correspond to SNM and PNM 
respectively. As expected, the binding energy per nucleon Eb of SNM initially decreases with increase in density, 
reaches a minimum a,t p — po and then increases. In case of PNM, the binding energy increases monotonically with 
increasing density in consistence with its well known behaviour. In Fig.[2Ka), we compare the Eb of SNM as a function 
of the nucleon density with a few representative results in the literature, namely, the Walecka model [l^] (long-short 
dashed curve), the DBHF calculations of Li et al. with Bonn A potential (short-dashed curve) (data for both the 
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models are taken from and the variational A18 + 6v + UIX* (corrected) model of Akmal at al. (APR) [2l| 

(long-dashed curve). While the Walecka and Bonn A models are relativistic, the variational model is nonrelativistic 
with relativistic effects and three body correlations introduced successively. Our model produces an EOS softer than 
that of Walecka and Bonn A, but stiffer than the variational calculation results of the Argonne group. It is well-known 
that the Walecka model yields a very high compressibility K. However, its improvised versions developed later with 
self- and cross-couplings of the meson fields have been able to bring down the compressibility modulus in the ball park 
of 230±10 MeV Q. Our model yields nuclear matter saturation properties correctly alongwith the compressibility of 
K = 260 MeV which is resonably close to the empirical data. In Fig. [D^b), we plot Eb as a function of the relative 
nucleon density for PNM. Similar to the SNM case, our EOS is softer than that of Walecka and Bonn A models, but 
stiffer than the variational model. We use this EOS to calculate the mass and radius of a neutron star of PNM as 
discussed later. 

The density dependence of pressure of SNM and PNM are calculated using the eqn. [551 These results are plotted 
(solid blue curves) in Figs.[3l^a) and (b). Recently, Danielewicz et al. have deduced the empirical bounds on the 
EOS in the density range of 2 < p/po < 4.6 by analysing the flow data of matter from the fireball of Au+Au 
heavy ion collision experiments both for SNM and PNM. These bounds are represented by the color-filled and shaded 
regions of the two figures. These bounds rule out both the "very stiff" and the "very soft" classes of EOSs produced, 
for example, by some variants of RMF calculations and Fermi motion of a pure neutron gas :4j. As shown in these 
figures, the EOS of SNM and PNM generated by our model are consistent with both the bounds. 




FIG. 2: (a)The binding energy per nucleon Eb as a function of relative nucleon density p/po for SNM. The results of present 
work (P.W.) are compared with the results of DBHF calculations with Bonn A potential [l^, the variational calculations of 
the Argonne group [Ij and the Walecka model '^2?]. The data for the Bonn A and Walecka model curves are taken from [l3| . 
(b) Same as Fig- (2a), but for PNM. 

The potentials per nucleon in our model can be defined from the meson dependent energy terms of eqs. (jl9l) . (j2ip 
and (|23p . Contribution to potential from the scalar part of the meson interaction is due to the pion condensates and 
is given by Vs = {hint + h^)/p as defined earlier. The contribution by vector mesons has two components, namely, 
due to the uj and the p mesons and is given hy Vy — + Vp = {h^ + hp)/ p. In the Figs. S] (a) and (b), we plot Vs 
and Vy as functions of relative density p/ po calculated for PNM (Fig. [D^a)) and for SNM (Fig.dJJb)) respectively. The 
magnitudes of the potentials calculated by our model are weaker compared to those produced by DBHF calculations 
with Bonn A interaction [l^ as shown in both the panels of Fig. H) In Fig. IDJa), we show the contributions to the 
repulsive vector potential due to uj mesons (short-dashed curve), p mesons (long-dashed curve) and their combined 
contribution (long-short-dashed curve). The contribution due to p mesons rises linearly at a slow rate and has a low 
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FIG. 3; (a) The pressure as a function of relative nucleon density for SNM as generated by the present work (P.W.) (sohd blue 
curve). The color-filled region in green corresponds to the bounds deduced from experimental fiow data and simulations studies 
by Danielewicz et al. I4j. The data for the curves corresponding to RMF(NL3) calculations and the variational calculations of 
Akmal et al. (APR) are taken from :4|.(b) Pressure as a function of relative nucleon density for PNM. The shaded region and the 
color-filled region in g reen correspond to the bounds deduced by Danielewicz et al. using the "stifi'" and "soft" parametrizations 
of Prakash et al. [33. Our EOS is consistent with these bounds in the cases of both SNM and PNM. 



contribution at saturation density. This indicates that major contribution to the short-range repulsion part of nuclear 
force is from oj meson interaction. 

Knowledge of density dependence of symmetry energy is expected to play a key role in understanding the structure 
and properties of neutron-rich nuclei and neutron stars at densities above and below the saturation density. Therefore 
this problem has been receiving considerable attention of late. Several theoretical and experimental investigations 
addressing this problem have been reported ([1, H, [s^ and references therein). While the results of independent 
studies show reasonable consistency at sub-saturation densities p < po, they are at wide variance with each other at 
supra-saturation densities p > pq. This wide variation has given rise to the so-called classification of "soft" and "stiff" 
dependence of symmetry energy on density jssl . [39| . 

Fig. [D shows a representation of the spectrum of such results alongwith the results of the present work (solid blue 
cu rve) . While the Gogny and Skyrme forces (dark rib-dotted and dotted curves respectively with data taken from 
[IjUll) produce "soft" dependence on one end, the NL3 force (dot-dashed curve with data taken from Q) produces a 
very "stiff" dependence on the other end. The analysis of experimental and simulation studies of intermediate energy 
heavy- ion reactions as reported by Shetty et al. [s^ (red triangles and long-short-dashed red curve repectively) , results 
of DBHF calculations of Li et al. and Huber et al. [H, H^, |40| (rib-dashed and magenta ribbed curve), variational 
model [1, [2l[ (short-dashed curve) , RMF calculations with nonlinear Walecka model including p mesons by Liu et 
al.^d^ (long-dashed green curve) as shown in Fig. [5] suggest "stiff" dependence with various degrees of stiffness. 
The experimental results (represented by the red triangles with data taken from Shetty et al. |39i|) are derived from 
the isoscaling parameter a which, in turn, is obtained from relative isotopic yields due to multifragmentation of 
excited nuclei produced by bombarding beams of ^^Fe and ^^Ni on ^^Fe and ^*Ni targets. Shetty et al. have shown 
that the results of multifragmentation simulation studies carried out with Antisymmetrized Molecular Dynamics 
(AMD) model using Gogny-AS interaction and Statistical Multifragmentation Model (SMM) are consistent with the 
above-mentioned experimental results and suggest (as shown by the red long-short-dashed curve) a moderately stiff 
dependence of the symmetry energy on density. Our results (represented by the solid blue curve) calculated using 
eqn. (P7|) are consistent with these results at subsaturation densities but are stiffer at supra-saturation densities. More 
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FIG. 4: (a) The potentials Vs, and Vp (as defined in the text) in PNM as calculated by our model are compared with 
the Bonn A results of Li et aZ. The contributions made by the tj-meson (short-dashed curve) and p- meson (long-dashed 
curve) mediated interactions are distinctly shown for comparison, (b) The potentials in SNM. Because of isospin symmetry, 
Vp (see text for definition) vanishes. Both the scalar (solid curve) and vector (short-dashed curve) potentials produced by our 
calculations are weaker in magnitude compared to those of Bonn A calculations. 



observational or experimental information is required to be built into our model to further constrain the symmetry 
energy at higher densities. In Figl^l the curve due to Huber et al. ^0] (with data taken from ^]) correspond to 
their DBHF 'HD' model calculations which involves only the a, lo and p mesons. Similarly the long-dashed green 
curve due to Liu et al. [s^ is from the basic non-linear Walecka model with a, uj and p mesons. Our formalism is the 
closest to these two models with the exception that in our model the effect of a mesons is simulated by the tt meson 
condensates. It is also noteworthy that our results are consistent with these results for densities upto 2pQ. 

The wide variation of density dependence of symmetry energy at supra-saturation densities has given rise to the 
need of constraining it. As discussed by Shetty et al [111, a general functional form E'gym = E^ymip/P^V 
emerged. Studies by various groups have produced the fits with E^y^^ '--^ 31 — 33 MeV and 7 ~ 0.55 — 1.05. A similar 
parametrization of the Esym produced by our EOS with E^y^ = 31 MeV yields the exponent parameter 7 = 0.85. 

We next use the equation of state for PNM derived by our model in the Tolman-Oppenheimer-Volkoff (TOV) 
equation to calculate the mass and radius of a PNM neutron star. The mass and radius of the star are found to be 
2.25M0 and 11.7 km respectively. 



IV. CONCLUSION 



In this work we have presented a quantum mechanical nonperturbative formalism to study cold asymmetric nuclear 
matter using a variational method. The system is assumed to be a collection of nucleons interacting via exchange of tt 
pairs, u) and p mesons. The equation of state (EOS) for different values of asymmetry parameter is derived from the 
dynamics of the interacting system in a self-consistent manner. This formalism yields results similar to those of the 
ah initio DBHF models, variational models and the RMF models without invoking the a mesons. The compressibility 
modulus and effective mass are found to be iiT = 260 MeV and M* /M — 0.81 respectively. The symmetry energy 
calculated from the EOS suggests a moderately "stiff" dependence at supra-saturation densities and corroborates the 
recent arguments of Shetty et al. j39| . A parametrization of the density dependence of symmetry energy of the form 
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FIG. 5: Symmetry energy Esym calculated from the EOS (as in Eq. I27|l (P.W.) (solid blue line) is plotted as a function of 
density along with results of other groups. The data for experimental points and the results of the antisymmetrized molecular 
dynamics (AMEy)_ simulations with Gogny-AS and Gogny interactions are taken from Shetty et al 39], DBHF (Bonn A) results 
are taken from 2^, RMF (NLp) data are from [30|] , the variational model of Akmal et al. (APR) 21] results are from ,3], DBHF 
{aujp) model of Huber et al. [40] data are from [29|, the Skyrme amd NL3 results are from ,8]. Our result shows consistency 
with those of other groups and corroborates the moderately "stiff" dependence of Esym as advocated by Shetty et al. Js^. 



Esym = E^ymiP/ PoV with the symmetry energy E^y^ at saturation density being 31 MeV produces 7 = 0.85. The 
EOS of pure neutron matter (PNM) derived by the formalism yields the mass and radius of a PNM neutron star to 
be 2.25Mq and 11.7 km respectively. 
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